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Abstract. This is a brief survey of certain constants associated with 
random lattice models, including self-avoiding walks, polyominoes, the Lenz- 
Ising model, monomers and dimers, ice models, hard squares and hexagons, 
and percolation models. 

1. Introduction 

Random lattice models give rise to combinatorial problems that are often easily-stated 
but intractable. This article briefly presents several such problems, each involving the 
numerical estimation of asymptotic growth constants. Emphasis is given to recent 
developments. The author has been collecting and writing about various mathemati- 
cal constants for five years. His large evolving website [18| provides further discussion 



and references; updates and corrections from readers are always welcome. 

2. Self- Avoiding Walks 
Let L denote the <i-dimensional cubic lattice whose sites (vertices) are precisely all 
integer points in <i-dimensional space. An n-step self-avoiding walk uj on L, be- 
ginning at the origin, is a sequence of sites u(0), uj(1), u(2), u{n) with u(0) = 0, 
\uj{j + 1) — u>(j)\ = 1 for all j and uj(i) ^ for all i ^ j. 

The number of such walks u is denoted by c(n). For example, c(0) = 1, c(l) = 2d 
and c(2) = 2d(2d- 1). 

What can be said about the asymptotics of c(n)l It is known that 



u = lim c(n) 



exists and is nonzero. The value fi is called the connective constant and clearly 
depends on d. The current best rigorous lower and upper bounds for /i, plus the best 
known estimates for /z, are given in the following table |35| , [|IJ, |37| , JR^, p3| : 



d 


lower bound 


best estimate 


upper bound 


2 


2.62002 


2.6381585 


2.6939 


3 


4.572140 


4.683907 


4.7476 


4 


6.742945 


6.7720 


6.8179 


5 


8.828529 


8.8386 


8.8602 


6 


10.874038 


10.8788 


10.8886 



1 
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The upper bounds were computed only recently via the Goulden- Jackson cluster 



method fl38fl . Similar techniques can be used to estimate other constants associated 
with the combinatorics of words, e.g., the asymptotics of ternary square-free words 
and of binary cube- free words [Q . 

The connective constant values fi given above apply not only to the growth of 
the number of self- avoiding walks, but also to the growth of numbers of self- avoiding 
polygons and of self-avoiding walks with prescribed endpoints [[2711 . 

For d = 2 and 3, there apparently is a positive constant 7 such that 

lim C(n) 



exists and is nonzero |33|], JT(J, fl2f . These constants are conjectured to be || if d = 2 



and 1.1575... if d = 3. 

Another interesting object of study is the mean square displacement 



S (n) = £{Kn)| 2 } = — -$>H| 



c{n) 

where the summation is over all n-step self-avoiding walks to on L, weighted with 
uniform probability. Like c(n), it's believed for d = 2 and 3 that there is a positive 
constant v such that 

lim S{U) 



exists and is nonzero; moreover |$3|, [|K|, it's conjectured that v — | if d — 2 and 
v = 0.5877... if d — 3. The critical exponents 7 and v are thought to be universal 
in the sense that they are lattice-independent (although dimension-dependent). No 
one has yet discovered a proof of their existence on L, let alone a proof of universality. 

3. POLYMINOES 
A domino is a pair of adjacent squares. 

<FIGURE 1> 

Generalizing, a polyomino or lattice animal of order n is a connected set of n 
adjacent squares, e.g., for n = 3, 

<FIGURE 2> 



Define A(n) to be the number of polyminoes of order n, where it's agreed that two 
polyminoes are distinct iff they have different shapes or different orientations. There 
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are different senses in which polyminoes are defined, e.g., free versus fixed, bond 
versus site, simply-connected versus not necessarily so, and others. For brevity's 
sake, we focus only on the fixed, site, possibly multiply-connected case. 



Redelmeier |43| computed A{n) up to n = 24 and Conway and Guttmann [16| 
extended the sequence to n = 25. It is known that the limit 

a = lim A(n)™ 

n— >oo 

exists and is nonzero. The best known bounds on a are 3.791 < a < 4.649551 
as discussed in [^9J, J4"2"| , [ f4"9"| , [ f2T| . Improvements are possible using the new value 



A(25). The best known estimate, obtained via series expansion analysis by differential 
approximants ||16||, is a = 4.06265... A more precise asymptotic expression [ITO] for 



A(n) is 

C ■ rT l ■ a n 

where C = 0.316..., but such an empirical result is far from being rigorously proved. 

4. Lenz-Ising Model 
Let Lis; denote the regular d- dimensional cubic lattice with N = n d sites. For exam- 
ple, in two dimensions, is the n x n square lattice with = n 2 . To eliminate 
boundary effects, Ln is wrapped around to form a <i-dimensional torus so that, with- 
out exception, every site has 2d nearest neighbors. 

Let's agree that a nonempty subgraph of is connected and contains at least 
one bond (edge). Suppose that several subgraphs are drawn on L N with the property 
that 

• each bond of is used at most once 

• each site of is used an even number of times (possibly zero). 
Call such a configuration on Ljv an even polygonal drawing. 

<FIGURE 3> 

Note that an even polygonal drawing is the union of simple, closed, bond-disjoint 
polygons but need not be connected. Other names in the literature for these con- 
figurations include closed or Eulerian subgraphs. Let B(r) be the number of even 
polygonal drawings for which there are exactly r bonds. For example, when d = 2, 
it follows that B(A) = N, B(Q) = 2N and B(8) = ±N(N + 9) . When d = 3, it 
follows that 5(4) = 3N, B(Q) = 22N and B(8) = \N(9N + 375) . Computing B(r) 
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for larger r is quite complicated, especially considering that in > 2 dimensions, the 
drawings can intertwine and be knotted! 

Define the high temperature zero magnetic field free energy for the Ising 
model to be the series 

1 CO 

iv r fc=o 

where z is called the activity. If T denotes temperature, then tanh - (z) oc 1/T . 
Coefficients (3 k of this series, as polynomials in d, are given here : 

k P k 

4 \d(d-l) 

6 |d(d-l)(8d-13) 

8 \d\d- l)(108rf 2 -424c/ + 425) 

10 ±d(d - 1) (2976d 3 - 19814c/ 2 + 44956ci - 34419) 



The radius of convergence 0, 0, |§, |20|, [|IJ for /3(z) is 

f V2-1 



lim /3 2 21 



k— too 



'2k 



0.414213562373095049... if d 

0.218094... if d 

0.14855... if d 

0.1134... if d 

0.0920... if d = Q 

0.0775... if d = 7 



2 
3 
4 
5 



which is important since knowledge of z c gives the critical temperature or Curie point 
T c of the model. The exact two dimensional result is a famous outcome of work by 
Kramers and Wannier and by Onsager. 

Here is a related problem. Suppose that several subgraphs are drawn on L N with 
the property that 

• each bond of Ln is used at most once 

• all sites of Ln, except two, are even 

• the two remaining sites are odd and must lie in the same (connected) subgraph. 

Call this configuration an odd polygonal drawing. Note that an odd polygonal 
drawing is the bond-disjoint union of an even polygonal drawing and an (undirected) 
self-avoiding walk linking the two odd sites. 



<FIGURE 4> 
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Enumerating odd polygonal drawings gives rise to what physicists call the high tem- 
perature zero field magnetic susceptibility 



z 



A- 



k=0 



Coefficients Xk °f this series, as polynomials in d, are listed in [p0| . The radius of 



convergence of x( z ) is the same as that for (5{z) for d > 1. Further, when d = 2, 



7 



lim (1-7^ ■ xO) = 0.9625817322. 



T 



Wu, McCoy, Tracy, Barouch |52] , (T^j determined an exact formula for this coefficient 
in terms of the Painleve III function. The expression is complicated: can a simpler 
formula in terms of other mathematical constants (e.g., Glaisher's constant [|18[]) be 
found? An exact analog of this formula for d — 3 is also evidently not known. 

5. Monomers and Dimers 
Let Ljy be as before, but without wraparound. Two sites of are called adjacent 
if the distance between them is 1. A dimer consists of two adjacent sites of and 
the (non-oriented) bond connecting them. A dimer arrangement is a collection of 
disjoint dimers on Ljv • Uncovered sites are called monomers, so dimer arrangements 
are also known as monomer-dimer coverings. A dimer covering is a dimer 
arrangement whose union contains all the sites of L^. 

For d = 2, let g(n) denote the number of distinct monomer-dimer coverings of 
L N , then clearly g(l) = 1, g{2) = 7 and asymptotically ||, |25|] 



K = lim g(n)N = 1.940215351... 

No exact expression for the constant k is known. Baxter's approach for estimating 
k was based on the corner transfer matrix variational approach. A natural way for 
physicists to discuss the monomer-dimer problem is to introduce an activity z for the 
number of monomers. The constant k then corresponds to the situation in which 
z = 1. Values of g(n) were recently computed p6] up to n — 21. 

The above contrasts with the special case of dimer coverings. An exact expression 
is known here for d = 2, due to Kastelyn, Fisher and Temperley. If f(n) is the number 
of distinct dimer coverings of L^, then f(n) = if n is odd and asymptotically 

2 2G 

Jim fin)* = exp(— ) = 1.79162281206959342... 



n even 
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a fascinating and unexpected occurrence of Catalan's constant G. For d — 3, the 



number h(n) of distinct dimer coverings [|32| of Ljy is h(2) = 9, h(4) = 5051532105 
and asymptotically 

2 

lim h(n)N = exp(A) 

n even 

where the constant A is known only imprecisely. The current best rigorous bounds 



HI, |J, @] are 0.44007584 < A < 0.463107 and the best known estimate § is 



A = 0.4466... . 

6. Ice Models 

Let Ljv be as before, with wraparound. An orientation of is an assignment of 
a direction (arrow) to each bond of L^. Assume that d = 2 henceforth. What is 
the number, 8(n), of orientations of such that at each site there are exactly two 
inward and two outward pointing edges? (Such orientations are said to obey the ice 
rule and are also called Eulerian orientations.) Here is a sample configuration: 



<FIGURE 5> 



After intricate analysis, Lieb proved that [ 40fl , |45 



lim 6(n) 



4X§ 



1.539600717839002039... 



This constant is known as the residual entropy for square ice. Interestingly, 36 (n) 
is also the number of ways of coloring the square faces of with three colors so that 
no two adjacent faces are colored alike [|J. Values of 6{n) were recently computed 
26| up to n = 13, and relevant computational complexity issues were discussed in 



|35 

The residual entropy W for ordinary hexagonal ice Ice-Ih and for cubic ice Ice-Ic 
(both complicated three-dimensional lattices) satisfy 



1.5067 < W < 1.5070 



and are equal within the limits of Nagle's estimation error. These configurations 
are not the same as the simple models mathematicians tend to focus on. It would 
be interesting to see the value of W for the customary n x n x n cubic lattice L3, 
either with the ice rule in effect (two arrows point out and four arrows point in) or 
with Eulerian orientation (three arrows point out and three arrows point in). No one 
appears to have done this. 
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7. Hard Squares and Hexagons 

Consider the set of all n x n binary matrices. What is the number F(n) of such 
matrices with no pairs of adjacent l's? Two l's are said to be adjacent if they lie in 
positions and (i + or if they lie in positions and + 1), for some i, 
j. Equivalently, F(n) is the number of configurations of non- attacking Princes on an 
n x n chessboard, where a "Prince" attacks the four adjacent, non-diagonal places. 



Let N 



n 



then 



£ = lim F(n 



N 



1.50304808247533226... 



is the hard square entropy constant ||, [0, j34), ||. Essentially nothing is 
known about the arithmetic character of £. 

Instead of an n x n binary matrix, consider an n x n binary array which looks 
like: 



/ a u 



021 



031 



041 



V 



022 



032 



O42 



052 



023 



033 



O43 



053 



\ 



O34 



Q.44 



O54 



064 / 



(here n — 4). What is the number G(n) of such arrays with no pairs of adjacent l's? 
Two l's here are said to be adjacent if they lie in positions and (i + 1, j), or in 
and (i, j + 1), or in and (i + l,j + 1), for some i, j. Equivalently, G{n) 
is the number of configurations of non-attacking Kings on an n x n chessboard with 
regular hexagonal cells. It's surprising that the hard hexagon entropy constant 



77 = lim G(u)n = 1.395485972479302735... 



is algebraic (in fact, is solvable in radicals |28| ) with minimal integer polynomial |5 

25937424601x 24 + 2013290651222784x 22 + 2505062311720673792x 20 + 
797726698866658379776x 18 + 7449488310131083100160x 16 + 



2958015038376958230528x 14 - 72405670285649161617408a: 1 - + 
107155448150443388043264x 10 - 71220809441400405884928x 8 
-73347491183630103871488x 6 + 97143135277377575190528a; 4 



,12 



-32751691810479015985152 



This is a consequence of Baxter's exact solution of the hard hexagon model , [0 via 
theta elliptic functions and the Roger s-Ramanuj an identities from number theory! 
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Just as series for the Ising model were defined using counts of even polygonal 
drawings with exactly r bonds, series for the hard hexagon model can be defined 
using counts of non-attacking configurations of exactly r Kings (and likewise for the 
hard square model). The radius of convergence for the hexagon series 

11 4- 5^/5 

z c = — = 11.09016994374947424... 

2 

possesses an exact expression. No similar theoretical breakthrough has occurred for 
the square model, hence the radius of convergence for the square series 

z c = 3.7962... 

has no known analogous formula ||. These values are important since they corre- 
spond to the critical activity (e.g., temperature or density) at which a phase transition 
occurs in the model. 

If one replaces Princes by Kings on the chessboard with square cells, then the 
corresponding constant p4]| is 1.342643951124.... A related problem, that of enumer- 



ating the maximal configurations of -j nonattacking Kings, was discussed in |)T] 







8. Percolation Models 
Let M be a random n x n binary matrix satisfying 

• rriij = 1 with probability p, with probability 1 — p for each i, j 

• rriij and mu are independent for all ^ (k,l). 

An s-cluster is an isolated grouping of s adjacent l's in M, where adjacency means 
horizontal or vertical neighbors (not diagonal). For example, the 4x4 matrix 



M 



has one 1-cluster, two 2-clusters and one 4-cluster. The total number of clusters is 
4 in this case. For arbitrary n, the total cluster count is a random variable with 
normalized expected value 

E(C n ) 








1 


M 


1 


1 











1 





1 


V i 








1 / 



n 2 



The limit K${p) of K n exists as n approaches infinity, and Ks(p) is called the mean 
cluster count per site or the mean cluster density for the site percolation 
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model. It's known that Ksip) is twice continuously different iable on [0,1]; further, 
Ksip) is analytic on [0,1] except possibly at one point p = p c . Monte Carlo simulation 
and numerical Pade approximants can be used to compute Ksip)- F° r example [54]. 
it's known that K s i\) = 0.065770.... 

Instead of an n x n binary matrix M, consider a binary array A of 2n(n — 1) 
entries which looks like 



A 



V 





«12 




ai4 




ai6 




an 


«22 


«13 


^24 


«15 


«26 


ax? 




«32 


^23 


«34 


^25 


a 36 


Ct27 


0-31 


«42 


^33 


O44 


^35 


«46 


^37 



/ 



(here n — 4). One should associate a^- not with a site of the n x n square lattice 
(as one does for m^) but with a bond. An s-cluster here is an isolated, connected 
subgraph of the graph of all bonds associated with l's. For example, the array 



.4 



V 



1 









1 
1 












has one 1-cluster, one 2-cluster and one 4-cluster. For bond percolation models 

such as this, one often includes 0-clusters in the total count as well, that is, isolated 
sites with no attached l's bonds. In this case there are seven 0-clusters, hence the 
total number of clusters C 4 is 10. The limiting mean cluster density 

E{C n ) 



K B ip) 



lim K n 

n—too 



lim 

n— >oo 



rr 



exists as n approaches infinity and similar smoothness properties hold. Remarkably, 
however, an exact integral expression exists at p = \ for the limiting mean cluster 
density M, W7 



K 



1 d 

r ot{y) -d-y 




sech 



7TX 



•In 



'cosh (x) — cos(2y) 
cosh (x) — 1 



dx 



V= a 
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which Adamchik [53] recently simplified to 

3\/3-5 



*4> 



0.09807621135331594... 



An analogous expression for the limiting variance of cluster density was computed in 



5j| . In the same way, the bond percolation model on the triangular lattice gives a 
known limiting mean cluster density at a specific value of p (discussed below), but 
the relevant variance is not known here. 

Let's turn attention away from mean cluster density K(p) and instead toward 
mean cluster size S(p). In the examples given earlier, S4 = (1 + 2 + 2 + 4) /4 = 9/4 
for the site case, S4 = (1 + 2 + 4)/3 = 7/3 for the bond case, and S(p) is the limiting 
value of E(S n ) as n approaches infinity. The critical probability or percolation 
threshold p c is defined to be 



Pc 



inf p 

0<P<1 

S(p) = OD 



that is, the concentration at which an infinite cluster appears in the infinite lattice. 



For site percolation on the square lattice, there are rigorous bounds | 5"0fl , [IE 

0.556 <p c < 0.679492 



and an estimate [53| p c = 0.5927460.... based on extensive simulation. Ziff [5_4| has 
additionally calculated that Ks(p c ) = 0.0275981... via simulation. 

In contrast, for bond percolation on the square and triangular lattices, there are 
exact results due to Sykes and Essam. Keston |27j proved that p c = ~ on the square 
lattice, corresponding to the expression Kb{\) above. On the triangular lattice, 



Wierman [27] proved that 



Pc 



, 7T 

2 ■ sin — 

,18 



0.347296355333860698. 



corresponding to another exact expression M, [54 



35 3 / vr \ 23 3 

K B (p c ) = -esc — = 

yF ' 4 2 U8/ 4 2 

= 0.1118442752845497... 



vY4 • (l + ^v^) + VV4 • (l - iV3 



Similar results apply for the hexagonal lattice by duality. 

Much energy has been placed into the computation of universal exponents for 
random lattice models in <i-dimensional space (akin to what we briefly mentioned for 
self-avoiding walks). The existence of such exponents is hypothesized on the basis of 
both theoretical procedures (finite size scaling and renormalization group methods) 
and experimental data. We shall not attempt to discuss this important subject, but 



instead refer readers to J33|, [0, |27 [. 
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